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Abstract 



We re-derive hydrodynamical equations in General Relativity (GR) in the comoving reference frame 
for spherical symmetry and obtain from them the well-known but not explicitly derived Lagrangean 
equations in Special Relativity (SR), that is, for the case of weak gravitational fields. Explicit 
formulae are presented which relate General Relativistic independent variables, Lagrangean mass 
coordinate and variables in the lab-frame. Conversion of one set of variables into another requires 
knowing the solution to the set of equations. These formulae allow one to translate the solution to 
the exact set in General Relativity into the form in which Special Relativistic solutions are usually 
obtained. This is applicable for comparison of SR-numerical simulations of collapses, GRBs and 
Supernovae explosions with more precise GR-simulations. 
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1 General relativistic hydrodynamics in the comoving ref- 
erence frame 



A detailed derivation of the hydrodynamical equations for the spherically symmetric motion in 
3+1 dimensions can be found in the paper [5] (Liebendorfer, Mezzacappa & Thielemann) , along 
with relativistic Boltzmann equation and its forms. However, the relation of General Relativistic 
hydrodynamical equations with Special Relativistic ones, such as used without derivation by Daigne 
& Mochkovitch [2], has remained unclear. We fill this blank with a comprehensive derivation of the 
basic equations in [2] starting from the very Einstein equations. 

We will work in the reference frame moving with (baryonic) matter and suppose that the metric 
tensor has been diagonalised, adopting the following (already spherically symmetric) form of the 
interval: 

ds 2 = e 2 ^dt 2 - e 2X dR 2 - r 2 dQ 2 , (1) 

where dQ 2 = d6 2 + sin 2 9d(p 2 is a solid angle differential. The quantities (p, A, r are the functions 
of independent variables R, t, which are correspondingly the arbitrarily chosen radial and universal 
time coordinates. We have already set c to be 1. 

Jacobian of these 4-coordinates is y^g = e v+x r 2 \ sin6>|, and jacobian of spatial coordinates 
equals ^fy = e A r 2 |sin#|, hence, a 3-dimensional infinitesimal volume element integrated over ip, 9 
equals 

dV = eHnr 2 dR, (2) 

and r(R,t) has the geometrical sense that the length of circle of radius r is 2nr. 

Now we have to put some physical sense to the coordinate R in order to fix its definition (one 
can apply a conversion R —> R staying in the comoving frame and only changing A — * A, so R is 
not fixed). We relate this coordinate with the baryon density, so that dR/G = dm b is the baryon 
mass inside the layer dR: 

pbdV = pbe x 4irr 2 dR = drrib = —dR. (3) 

In the following we put G — 1 (so that mass and length are of the same dimension) and denote by 
prime the d/dR derivative. 

In our hands are the baryon number conservation law 

(P^ M );m = 

and the Einstein equations 

In our comoving reference frame the 4- velocity equals w M = (e _v ,0), and using the formula for 

= ^ we get 

(PbU^ = (Pbu^ + r^p b u» = e-*(p + p\ + 2p r -) = 0, 



that is, 



(We omit the index b of pb and m b ). The Ricci tensor for the diagonal metrics 

9ik = eie 2Ft 5 ik , e { = (1, -1, -1, -1), 



can be easily calculated using formulae 



m^i,l 



Rik — X! \_ F i,kFk,i + F,kF,i — F,iF,k — F,i,k , i ^ k. 



Finally the Einstein equations yield: 



-8*7? = e- 2A L" + V a + j + ±(<f' - A')) - e- 2 - (X + A 2 + f - + *-(\ - 0)) 



8ttT/ = -e- 2X 2— + — - 2— I + e 



r'X 



-2<p 



r 2 Xr\ 1 



(5.1) 

(5.2) 

(5.3) 
(5.4) 



8nT t R = e~ 2X - (f - <p'r - Ar') . 
The hydrodynamical momentum-energy tensor 

Tvu, = (e + p)u^u v - pg^ 
has very simple mixed components in the comoving frame: 

7? = e, T* = -p,i = 1,2,3, Tj = 0,i^j 

which should be inserted in the LHS of eqns. (5) 

Let us also write down the conservation laws (which follow from Einstein eqns(5) because the 
4-divergence of Einstein tensor equals 0) 



7^ = e + (e+p)(A + 2-)=0, 

T^ = -p'-<p'(e + p) = 0. 
We used here the following useful formula for divergence of a tensor: 



(6.1) 
(6.2) 



v -g v - 9 
From the comparison of eqns(6.1) and (4) we see that 



e + p p 



(7) 



which is nothing more than the adiabaticity relation dE +pdV = 0. Using this relation we also get 
the formula for the internal energy e in defined as e = p(l + e in ): 



d e d 1 

e in = ^t(-) = -P^t(-)- 



Now let's define the following quantities: 

f i 9r 
rh= 4irr e^—dR, 
J on 



(8) 



(9) 



U = e-vf, (10) 

r = e -V. (ii) 

From this definition we simply find using the relation (6.1) for e with excluded A via (5.4) and 
expressing ip' through (6.2): 

in = —4irr 2 pr (12) 
Equation (5.3) with T\ = e provides further relation between m(R,t) and T, U. Really, 

87ier 2 r' = (r - e" 2 V 2 + e~ 2v rf 2 )', 

and hence 

R 

2m(R, t)= J Mr 2 r'edR = r(l + U 2 - T 2 ). (13) 
o 

We do not have to forget the relation (3) which reads 
or equivalent ly 

r = - = 47rr 2 e A . (3") 
Now we are able to formulate equations for evolution of p, U,Si n . Using (7) and (6.1), 

. _ d 1 1 (it r'-r<p' \ A 2 (2r r' - t<p' \ e x 2 . , 1 2 , 

r = - — = = e 47rr = — (Anr re v ) e v = — (47rr (7 e% 

^ p p V r r' ; V r r' J r' V ; r J ' 

i.e. finally 

= ^r 2 ^)'. (14) 



Going further, 



r '2 i -2 

e^f) = e-^r - r<p) = e - 2A ( — + 2r>') - 47rrp e -2 * — , 

2?~ 27 1 27 1 



where we have used (5.1) for expressing (f — r<p)e 2lp , and finally 



777- r 777- 

e-viJ = e~ 2X 2r'(p' - 4?rrp = --4vrrV ? - 47rrp, (15) 



where we used the relation (13) for m/r and introduced the specific enthalpy 

h=(e + p)/p. (16) 

The set of equations (14,15,8) can be used as hydrodynamical equations without ultrarelativistic 
particles' fluxes (i.e., additional terms in the energy- momentum tensor). Let me mention that 
Einstein equations (5.2,5.3) were not used directly in the derivation, but instead two conservation 
laws (6.1,6.2) were used which are themselves consequences of the whole set (5). 

Now let's add the momentum-energy tensor of ultrarelativistic particles, let us call them "neu- 
trino" (but these may also be photons) 

T$ = {pE u + p u )u,u u + {F^u v + F^u„) - p u g, v . (17) 

p u is the neutrino pressure (not a 4- vector). 



We may calculate the divergence of a slightly more general tensor 

Tfiv = Au^u u + (BpU,, + B v u^) + Cg^, (18) 

which yields 



7% u = A + C + (*+\)A + 



/ • 2r • \ 2r' 

[Bq + B (- + A) J (1 + e"*) - B^ + B' R + B R (— + tp' + A') 

(19.1) 



T£.„ = C" - ^(A + 2B ) + e-* + B R (y + A)) . (19.2) 

Here 5^ = PrrB r = -~f RR B R = -e 2X B R . In our case A = e+p+pE u +p u , B^ = Fj?\ C = ~{p+p v ). 
We also introduce notations Q, q according to: 

= (e*Q, -e x q, 0, 0), T%W = (e~*Q, e~\ 0, 0), (20) 

so that Q, q may be expressed through e v = pE v , p v , F u = ^rrF^ 11 (4- vector of flux must have 
F = in the local rest frame): 

e^Q = e u + (e„+p„)(\ + - F' v - F v (<p' + A' + ^), 



e\ = p' v + y\e v + Vv ) - e -f (p v + F u {\ + —)); A + — = --. 

V r J r p 

The components of momentum-energy tensor are now 

T* = e + pE„; T t R = e v F^ R = e^ 2X F u ; T R = T e 9 = -(p + p v ) . 
From (5.4) we obtain now 

■^r'-^- 47 rre^ 
r' 

and (6) are replaced by 

£ = -(e + p)(A + ^)-e v g, (6.1") 

ip* = - P ' + ^ = - T -p> - — q -. (6.2") 
Some minor changes in other relations: 

hp = e + e^Q (7") 

e^e in = -pe'^r - rQ (8") 

m(R,t) = J ' 4nr 2 Tfr'dR (9") 

Equations (10,11,13) are not modified, but (14,15) acquire additional neutrino couplings in the 
RHS: 

e-*T = ^(4nr 2 U)'-T^F u , (14") 

r T 7TI 

e- v U = --4irr 2 (p' + — —q) - — - Anr(p + p v ). (15") 



2 Lagrangean special relativistic hydrodynamics 



After having formulated hydrodynamical equations in the comoving RF, we need to relate some 
quantities with observable values in some fixed RF. We'll do this in the limit of very small gravita- 
tional mass, when a pseudo-euclidean (Minkowskian) RF exists: 



ds 2 = dt 2 - dr l - r z dn 2 . (21) 

It is obvious that r here is the same as in the comoving RF (1), because the length of a circle (27rr) 
is the same in both RF's (Lorentz transformation does not change transverse sizes, and matter 
moves only radially). For the same reason dfl 2 is the same in both RF's, because it is the area on 
the unit sphere which is not transformed being transverse to the motion. 

We need to relate times t and t, and this must be some function t(t,m). Derivatives of this 
function can be easily determined. As for a moving particle (which has to = const), its worldline 
interval equals ds = e v dt, and on the other hand, 



( 8r\ ~ ~ 1 

ds = , dt 2 - hH dt 2 = dtVT^ = -dt, 

\ \ dt )m 7 



where we've defined v = (J^J , which is simply radial velocity of a particle (to) in the fixed RF, 
and usual Lorentz factor 7 = ■ Hence, 



I) =7e^ } . (22) 



Now determine (jj?) . For a stationary point (r = const) we have ds 2 = dt 2 = e 2ip dt 2 — e 2X dm 2 , 

dr/dt c 
dr/dm 



dm = — Qr/fhn dt (derivative along the curve r(m,t) = const). Therefore, 




r-*$£f- (23) 

And with this formula it is straightforward to find (^) : as (g) = (§f) + (j£) t (^jr) , 



(s t ^(ei-(a)^(^-v--m> ™ 

We recall that U — 71* (this follows from the equality U = (^j , ds = ^dt for a moving particle). 
Moreover, in the case of weak gravitational fields which we assumed, 



r = Jl + U 2 - — ^7- — « 7 . (25) 
V r 7 r 

The difference between T and 7, that is, is always small in the ultrarelativistic case: it goes like 
~ Ptot r ' 2 /l for small r, equals G ^ ot outside of the most mass of matter, and never exceeds about 



~ 5ep. So, derivatives (23,24) take the form: 



And when one wishes to express (t, m)-derivatives through (t, m)-derivatives, the next formulae are 
to be applied: 



d 



d 



75T = ^ e 



dt 



dm 



dm 



dt 



+ e x U 



dt 



(27.1) 
(27.2) 



We also write down for reference the formula relating distances in the two RFs: 

1 



dr\f 



ry2 



dr\ t . 



(And remember that physical length between two mass points in the comoving RF is dL = e x dm = 
^dr\ t ; it is now clear that the distance between the same points measured at the moment t in the 
euclidean RF is just dL/j.) One is able to see now, that from the m-derivatives of pressure come 
also t-derivatives, and this circumstance complicates the set of equations. For example, rewrite the 
equation (14) for density evolution (for simplicity without photon/neutrino flux) in (t, m)-terms: 



dr 4n 
1 dt = ~ 



4tt 

7 



(28) 



We can express here ^ using the second equation of our set (15), which reads in new variables 



dU 

7— T = -Anr 
dt 



p + -r 
h 



dm 



dt 



(29) 



We see that some additional t-derivatives appear on the RHS of evolutionary equations for (inverse) 
density, but fortunab 
g?7 = vd^v), we find 



density, but fortunately this can be avoided writing equations for ^ instead of K Really, since 



d 1 



1 d 



dt 7p 7 2 dm 



(4:7ir 2 v ) = 



Aitr v d'j 2v 



+ 



7 3 dm p'-fr ' 



i.e., introducing V rn = and denoting d/dt and d/dm by dot and prime, respectively, 



T 



T 



fryr 



Substituting here 7' = rfvv ' and bearing in mind that 

(dr_\ = l_ eX = _J_ 
\dm)r 7 47rr 2 7p 

(but compare this with eqn. (11) — this is different derivative because it is taken at constant t) we 
finally see 

V m - An(r 2 vY = 0. (30) 

The same trick may be applied to evade d/dt— derivatives in the RHS of eqn. for momentum density 
U if one uses S m = hU instead of U : 



S rn + 4titV = 



iiirph 

7 



For reference let me write down also the very equation for U : 

U = 



4:iir 2 , . U 4irrp 



h 



Anr 2 p 



7 



(31) 



(29') 



Using E m = fry — - R as the third variable (instead of e in ) we also do get a total cancellation of 



7 p 

t— derivatives (one should also express p through U and 7 = vU): 

E m + 4ir(r 2 -fvp)' = -4tt^^. (32) 

Let me also mention that if we had not put G — 1, the only additional multiple in the set (30-32) 
would be G in the RHS of (31) and (32). This means that these additional (with respect to basic 
eqns. [2]) terms are of the order 0(^-) and therefore must be neglected since we have already 
omitted such terms, for example, in (25). Moreover, to take into account these corrections, coming 
from General relativity, one should change the very assumption of existence of the metrics (21). 

Finally, let's introduce one "new" independent variable, following the notations of the work [2]: 

dm = 7p A-Kr 2 dr\i = — Aiir 2 dr\ t . (33) 

7 

Clearly, dm\i = dm\ t , so that m is a baryon rest-mass Lagrangean coordinate. Actually, m(t,r) 
defined as 



m(t,r) = J 7p Anr 2 dr\i 



ro 



is the same baryon rest mass like in the comoving RF but expressed as a function of (t,r): 

m(t, r) = m(t, r). 

This implies that our notations are identical to those in [2]. From this definition of m it is easy to 
derive an equation for V m : 

( dm\ Anr 2 

d , 9 dr . , _ . dr 9 d 2 r . . . dr 9 dv . d , , 9 , 
Vm = r 2 — = 4vr 2rr— + r 2 ^—) = An(2rv— + r 2 —) = —(4nr 2 v). 

at dm dm atom dm dm dm 

3 Conclusions 

Let us summarize our results. If one is to consider General Relativistic effects, the exact set of 
hydrodynamical equations (14,15,8) can be used. In the weak field limit, one can apply the easier 
(and more suitable for comparison with observational data) set of Special Relativistic equations 
(30-32) which neglect 0(GM/R) effects but are formulated in the fixed pseudoeuclidean spacetime, 
call this "lab frame". But if an exact set is used, it is desirable nevertheless to present results 
in (r, t)-description, where t is the universal time in lab-(observer-)frame. To establish this, the 
functions t(t,m) and r(t,m) should be found, e.g., via integrating together with a set of hydro 
equations the following 4 equations: 

9i „ di Xtt 

7T+ =e ^i 7T = eU i 
dt dm 

dr _ dr _ ^ A 

dt ' dm ^ 

This Jacobian matrix has an evident remarkable symmetry. 

So, if one knows a GR-solution, for instance, a path of a shock wave, msw(t), then in the 
lab-frame the solution is a curve r = r(t,m sw (t)), t = t(t,m sw (t)) parameterized by t. 
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